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Cyclotron Resonance in Parabolic Quantum Wells in Tilted Magnetic Field 


T: Introduction 


Experimentally, the application of the additional parallel 
component of the magnetic field to the sample provides insight 
into very nontrivial behaviour of the  quasi-two-dimensional 
(Q2D) system and enables extraction of additional information 
on its properties. The method has successfully been used to 
study the combined and coupled intersubband cyclotron reso- 
nances Il in Q2D system, where the coupling of modes is expec- 
ted to occur. It was shown to give the new insight into the 
spin structure of states in the fractional quantum Hall effect 
[2-4] customarily treated as purely two-dimensional phenomenon 
in which case the parallel field: should have no effect. 
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The material engineering has reached now the state of art 
when the quantum wells can be produced with a given shape of 
confining potential. In a recent work [5] the quantum wells 
with parabolic confining potential V(z) œ z^ in the growth 
direction (by varying quadratically the composition x in 
Al Ga, As material from 0 for z = 0 to 0.19 at the well edges) 
have been studied. The precision of material processing when 
combined .with very precise optical experiments gives very 
detailed information the handling of which requires better 
theoretical treatment not only of the many-body but also 
one-particle properties. 

Recently the problem of a single electron moving in a 
quantum well confined by the parabolic potential V(z) and with 
magnetic field B applied at angles @ #0 with respect to the z 
direction has been solved for isotropic material analytically 
{6,7]. The exact solution gives the possibility of correct 
interpretation of the data outside perturbative regimes of 
small angles 9 or weak fields B. ; 

In particular, the origin of a gap in optical spectra of 
quantum well studied in tilted magnetic field has: been 
clarified [8]. The line shape, however, seen in (Di and also in 
{9] remains still unexplained. The satellite structure observed 
on both lines in the coupled resonance in tilted field 
configuration remains visible at @= 0 and may be due to spin 
dependent many-body interactions mixing the spin states in a 
given energy level resulting from hybridization of subband and 
Landau levels. In this work we shall neglect all the. possible 
many-body complications and will discuss the one-particle 
property of the parabolic quantum well in tilted magnetic 
field. Special attention will be focused - on the anisotropy of 
the system, and tne polarization effects. 

The organization of the paper is as follows. In Section 2 
we shall present our anisotropic model Hamiltonian and Sketch 
the solution of corresponding  eigenproblem. Some comments on 
the degeneracy spectrum and the effects connected with 
anisotropy of the mass tensor wil be made. The calculation of 


matrix elements entering the formula for oscillator strengths 
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is presented in Section 3. This is the main part of our work. 
The discussion of results as well as the proposals of new 
experiments and limited comparison with available experimental 
data is given in last Section. 


2. Model and (exact) solution 
A single electron Hamiltonian in a confining potential 
V(z) and magnetic field B described by vector potential Å in 


effective mass approximation is given by 


H = 2 S am, + Vz), (1) 


where m, are effective masses in three principal directions, by 


assumption coinciding with (x,y,z) coordinates. P Zi (p Kä 
eA), p, 7 -ih 2 V(z) = m O'z^/2, where © has been taken [7] 
i , 


in the form AQ = E (0 + 1 intersubband energy). We take B = 
(0,B.,B ) = B (0, sin@, cos6) and represent it by the potential 


Á in a Landau-tike gauge Ñ = B(zsin@- y cos@ 0, 0). This 
gauge is more convenient than that used in [7], as we can work 
only in the f representation without need of mixed 
representation. Tt leads very easily to the explicit 


expressions for the wave function also in r' representation. 

The resulting Hamiltonian commutes with the momentum 
operator P. Thus Pt hk. is the constant of motion, where Di 
is the wave vector of the plane wave exp(ik x). Taking all this 
into account and shifting the y coordinate to y =y —- 
1k /cos8, with Y = h/eB we get 


22 
2 2 L2 2 maz 
H, = m (eBz singe + eBy' cos&)*- a e- em A EE = 7 
x x “2m éy* 2m 2z 
y 2 
(2) 


This Hamiltonian represents two coupied harmonic oscillators. 


; 2 
The coupling term em zy‘sin(26) can easily be eliminated 
x 
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` by proper rotation of scaled coordinates vn Z and Yn y* by an 


angle ó, i.e. 


=- x 
y" = Ym y'cosó + Ym zsinó (3) 
na Ga y'sinó + Ym zcosó 


We get 


o rsin2e 
TAS m ———— MA (4) 
Gr — à (cose - eine 


where y? = m/m. and o = eB/Ym m. is the parameter which for 
isotropic system reduces to the cyciotron frequency 
€. = eB/m. a 

The spectrum of the problem consists of two harmonic 
oscillator branches z 


Be GN) + 1/2) he, ae (Qui os 1/2)hwo, P. (5) 


NM 


with frequencies 


w= (a a^» sin'e)cos'ó E G^cos^esin^ó d : G^ rsin26sin26 , (6a) 


D 


2 2 d. 2. ad 2 — 2 el 
w= (a + o sin 8)cos ó + e cos Osin ó — 


Mim 


&-ysin26sin26 d (6b) ` 


The eigenfunction of (1).can be written as 


i ik x 
Pe = INM k > = - a e $ (y"/A) $ (z'/A ), (7) 
$ x YL Di a Wi 


where ee = b/w, Li is the linear dimension of the system 


SCH 
in the x-direction and for periodic in x direction boundary 
conditions we have ky = en man = 0 42... Energies (5) are 
x 


degenerated, they do not depend on k „While the wave functions 
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do. To calculate the degeneracy let us note, that the position 
of the center (y 9 of both oscillators depends on the k e and 2 


Explicitly 
n 2 
ose : ky ztgó 
ot cosa Y 
(8) 
2 
e 3 Ze i Ka = zctgo 
o2 cos8 Yr 2 


Applying the proper boundary conditions (the center of 
oscillator should not be outside the sample) to (3) we get the 
degeneracy g of levels to be 


e Xx = BE F 
g h eBcos@ r eB. > (9) 


independently of the factor » The perpendicular component of 
the magnetic field controls the degeneracy of levels. 

Inspection of the results for isotropic medium m;-m in 
the formulae (4) and (6) shows that the ‘effect of anisotropy 
can be accounted for by replacing es eB/m and @ in isotropic 
expressions through 


& — o = o (y sin + cos?e)"* , (10a) 

c c e 

8 — 6 = arctg(ytge) e (10b) 
The cyclotron frequency in (10a) has a simple interpretation. It 


is just the cyclotron frequency of the bulk material with 
cyclotron effective mass m (8), given by 


m 
= € SE se [ zY sin*e + coste ]'* (11) 
D ym My LJ 


The dependence of Ə on GO for the anisotropy coefficients 
corresponding to PbSe, Si, PbTe and Ge shows Pigu. In 
calculations we have assumed that principal axis of the valley 
is along the z direction. 
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Formula (10b) shows that anisotropy effectively changes 
the tilt angle from 8 to larger values @ for kal and t^ 


smaller values for y < 1. 


Figure 1.The effective tilt 
angle e versus the actual @ 
for various ratio K = m/m, 
corresponding to known 
semiconductors. Curve 1 
is for Ge (K=19.5); 2-PbTe 
(K=10.0); 3-Si (K=4.8); 4- 
PbSe (K=1.6). 


E? “ 60 80 | 
8 [degrees] | 


The eigenenergies (5) for few first quantym numbers  (N,M) 
are shown in figure 2 as a function of applied magnetic field 
(w "eB/m) for tilt angle e = 45? (full lines) together with two 
dimensional calculations (dashed lines). As already mentioned, 
in a strictly two dimensional case the parallel field has no 
effect on the behaviour of the system and the lines cross. It 
has to be noted that the two-dimensional cyclotron frequency 
differs from the three dimensional e used here. In two 
dimensions only the perpendicular component of the B field 
enters the expression for the cyclotron energy. In this figure, 
as well as in all the forthcoming results of numerical 
calculations we have set m ^m -m.. 

Comparison of the above exact results with perturbative 
ones [8,10] can easily be done. It reveals that some of the 
effects obtained here, e.g. gap opening, could not be foreseen 


in perturbation theory. These will be discussed in Section 4. 
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tin Ewo 


Figure 2. Energies (5) are plotted against cyclotron 
frequency (full curves). The numbers refer to (N.M).The dashed 
curves represent the spectrum of the two dimensional sample. 


3. Cyclotron resonance 


The knowledge of spectrum and wave functions enables us 
the calculation of optical absorption spectrum. For light 
polarized in Xs direction the strength of the transition is 
given by the real part of the dynamical conductivity el 9), In 
linear response theory within dipole approximation it is given 
by [10] 


3 


2 
rip) a B.E -NO) e el erf 
beaatat * Se A, RA La Sib E - Pai, (12) 


with 
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wit EUM Lem dv. Iw zf (13) 
ue er Ey UND Ou NA 


where Ne is the electron concentration in oth level and v= 


is the velacity operator. 

The calculation of <w,]v,|¥,> is rather simple. To this 
end we shall transform the components of the velocity operator 
to the coordinate system denoted previously by  (x'.y'",z") (see 
equation (3)) and calculate the matrix elements entering (13) 
with help of wave functions (7). The general result (for  »*4 


is 
et 
EHM WEEN FAS sen Pye sinósine - cosécosé) 
x x m 3 
+ Ê u Z? pyp Sin8cosó + cosésiné) } à (14a) 
ART e ES ; . 
x N' M' |v. | k,.N, ND e Py uu: C06 ana teu tins. (14b) 
y " 
Mere: m | ; 
€x N' M* | v, Jk HMG H enn Py re POS ang cons). (14c) 
z 


where for ES y" 2",n = M,N, n' = Mi NI we have 


Kr AA Eeler |n EE + msg k asa) 


n'n-í 2 net 


- FOF su] (155) 


Ws. ud. [/2 
D'Aen ren LA |n» ih ( 5 CR 
Note again that in order to get matrix elements for anisotropic 
System from that for isotropic one we have to multiply 
"isotropic" formulae for <v> by corresponding Yn/m. and 
replace o and 6 by the values given on the right hand side of 
(10). 
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4. Results and discussion 


Having determined the energy spectrum and matrix elements 
we are now in position to calculate the absorption spectrum 
described by the real part of the conductivity tensor (12). 

The selection rules for the intra-band transitions may be 
obtained from (14). This equation shows that in the case of 
normal magnetic field (9 = 0} the light polarized in the x-y 


plane induces only the cyciotron transitions 
IN.M > + |NM + 1> while the light polarized in the 2 
direction couples to the intersubband transitions 


TATE > |N + LM». Above simple picture breaks down when 6 = D 
Due to the subband-Landau level coupling induced by n y the 
light polarized in the x-y plane .(the z-direction) may also 
excité the intersubband (cyclotron) transitions. Consequentiy, 
the cyclotron (intersubband) resonance spectrum splits into two 


D 

^ n ES Ex Jun d 1 is pic 

branches w C Nea m oo P, In the case of isotropic 
Mot 


effective mass Eqs.5 and 6 give 
D = > 
(wy = {2 + Yos pL SI + aute?sin®e] "ko (16) 


The dispersion of these two branches caiculated for ©@ = 30° 
45° and 60? shows Fig.3a. From this figure and Eq.l6 we find 
that behaviour of w and w is not consistent with perturbation 
theory for large angles © in following points: i) as already 
noticed in (7) with increasing ^ the lower branch (w ) tends 
to Q cos@ (and not to N as would simple perturbation indicate) 
while the upper branch (w°) tends to €. (not to 9 cose J. ii) 
the gap energy Æ = Kiw-w) takes on a minimum value equal 
Eo Sine when pom D cose (not when o Cose = (3, iil with 
increasing © the gap center E = bb + w)/2 lwsteosg Shifts 
from E,, to lower energy as (1*3cos*& C and not like cosé 
given by perturbation theory. 

Assume for simplicity that only ground level p.0 > is 
occupied by electrons. Then from (12-15) we find that in the 
case of the light polarized in the xc direction the intensity 
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of the low (high) frequency resonance is proportional to the 


oscillator strength fd Gne io: Substituting (14) into 


(13) we obtain, for the isotropic case, the following result: 


mat Gut) cosi - e), (17a) 
rM (ww)? sin (6 - 6) , (17b) 
foo. 7 COS 6, (17c) 
aes. = sin's ^ (17d) 
Ka 5 M. ? (17e) 


where the angie ó is defined by (4). 

We have caiculated the oscillator strengths (for light 
polarized in the x and y directions) as a function of 9. 
Results are presented in Fig.3b,c. We see from this figure that 
Similarly as in the perturbed theory the intensity of the 
branch o (w ) decreases (increases) as o. increases. However, 
the dependence of the relative intensity of the low and high 
frequency transitions on e. for given 8 reveals features that 
are not apparent in the results based on the perturbation 
theory. First of all we observe that even for normally incident 
light the shape of the absorption spectrum strongly depends on 
the light polarization. When the radiation is polarized in 
the x- direction the oscillator strengths of the two resonances 
are equal at Ox OQ. In the case of the y polarization the 


resonances have equal intensity at w= /(cos28)*? 3 6 < 45°. 


But when above condition is not fulfilled then f ES 
00,01 00,10 

Thus we see that the numerical value for the crossing point 
3 i E Ma DA : 

frequency (o ) determined by condition DT LE is 


independent on © only in the case of the x polarization. It is 
interesting to note that the crossing point (o ` ) determined 
by the condition of minimum separation between vw and w 
is equal to © cos@. From above considerations we may conclude 
that in general ^ < Ons o and the assumption commonly used 


in the interpretation of the experimental results, that 
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Figure 3. (a) Dispersion relation  (16),(b) oscilator strength 


Ën " (y? 
E HR (c) oscillator Steen D Ronm calculated for the 
different tilt angles (------ jo aff qe) gm. Lë y eee? rand 


(—— y. aereo" 


228 W. Okulski, KL Wysokinski.. . 


* 
w Ve 
E(u = Ww) = Elo - w )= E, is in principle correct only when 


6 is small. This conclusion is true also for nonpolarized 
light. In the case of the normal incidence the intensity of the 
|0.0» » JN.M> transition is then proportional to 


7 _ fœ ty? 
Poo; v KS E e (18) 
The dependence of f oos and e on w for different tilt 


angles is presented in the  Fig.4. 


0 WE 10 Ss 20 


25 3.0 
bé Eg 
Figure 4, The averaged oscillator strength (18) vs. o cal- 
culated for the same values of the tilt angle as in Fig.3. 


We see that, as in the case of the y-polarization the 
difference between ai and ^Q is positive and increases with 
increasing 6. The insert to this figure shows the height of 
experimental lines read off from Figure. 3 in [5] showing the 
far-inftrared magneto transmission measurements on  quasi-three 
dimensional modulation doped Ai-Ga-As structure. The asymmetry 
in the heights of the lines with respect to the crossing point 
is in agreement with our results. The direct comparison of this 
experiment with our theory is not legitimate as we have omitted 
very important aspects of the phenomena - the depolarization 


and the excitcnic effects which are known to strongly change 
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the frequencies of transitions. In perturbative regime it has 
been proved [10] that these effects merely lead to the 
renormalisation of frequencies without influence on the  B-field 
dependence. Preliminary calculations of the depolarization 
effect in nonperturbative regime seem to indicate more complex 
influence of the many body interactions on the absorption 
spectrum. The problem is being studied now. 

In summary, we have calculated the energy spectrum and 
velocity matrix elements for the parabolic quantum wells (with 
anisotropic effective mass tensor) placed in strong tilted 
magnetic field. We have found the pronounced polarization 
effects (see Fig.3) not predicted in perturbative approach. In 
view of this it would be very desirable to have experimental 
data performed with linearly polarized light ( in the x and y 
directions ) and at normal incidence. 
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